Specific Heat of the Quantum Bragg Glass 



cn 
o 
o 

oo 



> 
o 
o 
m 

O 



o 



Gregory Schehr,^ Thierry Giamarchi,'^ and Pierre Le Doussal^ 

'LPTENS CNRS UMR 8549 24, Rue Lhomond 75231 Pans Cedex 05, France 
^DPMC, University of Geneva, 24 Quai Ernest- Ansermet, CH-1211 Geneva, Switzerland 

(Dated: February 1, 2008) 

We study the thermodynamics of the vibrational modes of a lattice pinned by impurity disorder 
in the absence of topological defects (Bragg glass phase). Using a replica variational method we 
compute the specific heat Cv in the quantum regime and find Cv oc at low temperatures in 
dimension three and two. The prefactor is controlled by the pinning length. The non trivial can- 
cellation of the linear term in Cv arises from the so-called marginality condition and has important 
consequences for other mean field models. 

PACS numbers: 



Understanding the behavior of the specific heat in dis- 
ordered and glassy systems is a longstanding theoreti- 
cal challenge flj. Its low temperature dependence has 
been observed to be linear in a variety of glasses, in- 
cluding amorphous solids, disordered crystals and spin 
glasses 1^ 1^ 1^ . A phenomenological interpretation was 
proposed assuming the existence of two level systems ^ 
leading to such a linear behavior. Despite its remark- 
able success for many systems the range of applicabil- 
ity and microscopic justification of such arguments are 
still open questions. Analytical calculations based on 
models of quantum solids with structural disorder ac- 
count for higher frequency features such as the "boson 
peak" |^i7„.Sj but produce only phonon- like specific heat 
Cy ~ T^, except in the classical limit On the 

other hand new possibility arises from recent develop- 
ments of the mean field methods |llL to quantum 
spin glasses 0, . The glass phase is there described 
by a replica symmetry broken saddle point solution and 
evidence was found |l,'lll4| for a linear specific heat both 
in the p-spin and Heisenberg spin glasses. Since even in 
these solvable models the specific heat remains delicate 
to extract analytically, a later numerical study claims in- 
stead ~ T'^ and compares to experiments. It is 
thus still an outstanding problem to understand the be- 
havior of the specific heat already at the level of the mean 
field description of quantum glasses. 

In addition to structural and spin glasses, there has 
been considerable recent interest in pinned elastic sys- 
tems where disorder originates from substrate impuri- 
ties. Such systems cover a wide range of experimental 
situations both in the classical and quantum li mit, such 
as vortex lattices in superconductors 16, 0, Il9| . 
electron crystals [i^, |^ , charge and spin density waves 
[2^ . disordered liquid crystals [3, etc.. All these sys- 
tems are characterized by a competition between disor- 
der and elasticity, which leads to pinning and glassy be- 
havior. Although the question of positional order and 
correlations is now better understood, with reasonable 
agreement between theory and experiments, the behav- 
ior of the specific heat is still largely not understood 



and, ill some cases, affected by non equilibration effects 

El m m m m la, . 

It is thus crucial to develop a first principle method to 
compute the temperature dependence of the specific heat 
in such disordered elastic systems. Here we address this 
question for a Bragg glass system allowing both for ther- 
mal and quantum fluctuations ^3 ^"^^ obtain a quite 
general formula to compute the specific heat. We find 
that the term proportional to T in the specific heat van- 
ishes. This leads to a behavior of Cy in all dimensions 
between 2 and 4. The Larkin Ovchinnikov pinning length 
is found to control the coefficient. The result holds for 
a periodic object, i.e. a Bragg glass, and for interfaces 
with continuous degrees of freedom. We elucidate the 
general mechanism leading, in the mean field approach, 
to the cancellation of the term linear in T. We discuss 
how this property allows to extract such a term in re- 
lated quantum spin glass models. In a companion paper 
we investigate the classical limit and the applications 
to vortex lattices. 

Let us consider a collection of interacting quantum ob- 
jects of mass /i whose position variables are denoted by 
a TO-component vector Ua{Ri,T). The equilibrium posi- 
tions Ri in the absence of any fluctuations form a per- 
fect lattice of spacing a. Interactions result in an elas- 
tic tensor $0,^(9) which describes the energy associated 
to small displacements, the phonon degrees of freedom. 
Impurity disorder is modeled by a r (imaginary time) 
independent gaussian random potential U{x) interacting 
with the local density p{x) = — Ri — u(i?i,r)). 

We will describe systems in the weak disorder regime 
a/Ra ^ 1 where Ra is the translational correlation 
length, e.g. in a Bragg glass phase where the condition 
\ua{Ri,T) — Ua{Ri + a, t)| <C a holds, no dislocations 
being present. At equilibrium the system is described 
by the partition function Z = Tre"^^ = / DuDIie~^/^ 
with Hamiltonian H = iJph + H^is'- 



Hph — - 



1 rn{qy 



'^Ua{q)<^a,f3{q)uf3{~q) 
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Hdis= / d'^xU{x)p{x,u{x)) 



(1) 

and its associated Euclidean quantum action — 5'[w,n] = 
/^Xjin„(g)5rUa(g) - H. Here = /^^ denotes 
integration on the first Brillouin zone and all integrals 
over the imaginary time variable r go from to (3h, (3 = 
1/T being the inverse temperature. 

For simplicity we illustrate the calculation on a 
isotropic system with <^a,f3{<l) — cq^^ap and denote dis- 
order correlations U{x)U{x') — A{x — x'). The disor- 
der average is performed using the replica trick Z'' — 
J Dite"'^""/'' and integrating over 11, after some manip- 
ulations 1^ one obtains the following replicated action 
SeS = Sph + Sdis with: 

^ph = J d^xdT^ Y.^V.Uaf + ^{drUa? 
a 

Sd\s ^ y d'^xdrdr' ^ R{ua{x, t) - Uh{x, r')) 

ab 

R{u) = pI^^k cos{K -u) (2) 



where f2 is the volume of the system, V{B) = 
— Pq Aif exp(— i?i4r^/2) and a;„ = 2nn/f3h are the 
Matsubara frequencies. We have implicitly taken the 
limit fc ^ in : we denote G{q,ujn) = Gaaiq,^n) and 
parametrize Ga^b{q,^n) by G{q,u), where < u < 1, 
which is ojn independent. One has 



G{q,UJn 



c(g2+c^2/„2) + s + /(a;„) 



(6) 



-On,0 T- 

cq^ Uc cq'^ 



dv [o-]{v) 
cq^ + [<^]{v) ' 



Similarly we take hBab{T) = {[ua{x,T) — Uf,(a;, 0)]^)/m 
with B{t) and B{u) which is r independent. Uc and 
[a]{v) are defined in Ref. [ll jH)/{mn) is then 

calculated using the solution of the variational equa- 
tions [m . Isif where the best trial Gaussian action Q 
is obtained by breaking the replica-symmetry (RSB). Al- 
though the RSB scheme, and the behavior of B{u), de- 
pend in general on d and on V{B), Cy{T) depends on T 
only through Mc, S and B — B{u > u^)'- 



K 



Here v = \fcf\x is the pure phonon velocity and A^^ = 
/ d'^xe'^^'^ A{x) denote the harmonics of the disorder cor- 
relator at the reciprocal lattice vectors K, and pa ~ 
the average areal density. 

This quantum model is then studied using the Gaus- 
sian Variational Method (GVM) on the imaginary time 
action [ill Isij , implemented by choosing a trial ac- 
tion: 



-W"{B) 



E = cRZ 



(7) 
(8) 



where Rc is the Larkin length, and (0 is the so called 
marginality condition (MC) which automatically holds 
here for d > 2. We recall here the variational equation 
[ll| , relevant for the computation of {T) : 



So = 7^ I iZY. ^ab^aiq, ^n)ub{-q, -i^n) (3) together with O, where 
^1-^"' •'inb=i - 



where fc ^ is the number of replicas. 5*0 minimizes 
the variational free energy Fvar = Fq + -^{ScS - Sq)so, 
where Fq denotes the free energy calculated with the trial 
action 5*0. The specific heat is defined as follows 



/K) = T / dT{l-COs{L0,,T)){V'{B{T))-V'{B)) (9) 



(10) 



P^Jq cq^ + CLul/v^ + S + /(CJ„) 
5^ ^ Sr^ f 1 ^ cos(cj„t) 



P^Jn cq^ + civl/v^ + S + /(w„) 



a(r) 



-T 



d'^F d(H) 



QT2 



dT 



(4) 



where (...) denotes a quantum and thermal average and 
— disorder averaged. Introducing replicas, one finds 
W)/^ = UEaHphiua^Ila) + ^Sdis)s^„- Integrating 
over the Ha fields we compute the resulting average us- 
ing Sq (which amounts exactly - thanks to the variational 
equations - to compute Cu using the variational free en- 
ergy instead of the exact one). We obtain: 



1 c 

2 + 2 
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{q^ jW^)G'(g,iw„) 



The solution of these coupled equations can be organized 
as an expansion in powers of Ti, keeping pti fixed: S = 

0{n^) where one finds that the lowest order solution 
does not depend explicitly on (3% (of course I{ujn) al- 
ways depends implicitly on ph through ujn)- In general, 
{H)/n and Cy{T) will be functions both of h and (3h : 
(H)/{Tnn){h,ph) = Ho + nHi{ph) + h^H2il3h)..., where 
Hq is independent of Ph from which it follows that 

C,(T) = G,{h, (3h) = GoiPh) + hGiiph) + ... (11) 

After some manipulations, we find for d > 2 



dT{V{hB{T)) - / duV{nB{u))) 



(5) 



TiHi 



pn 



cq^ + So + lojuJn) 

q Cq^ + + So + Io{uJn) 



(12) 
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where in that Umit Sq is fixed by (O setting B = and 
/o(wn) satisfies the equation [lH 

/oK) - -4l^"(0)(Ji(Eo) - Ji(So + + ^oK))) 



g (Cg2 + ZY 



(13) 



from which we can extract the low lo behavior of its ana- 
lytic continuation Io{LUn —iuj + O'^) = ^o(^) +*^o'('-^) • 



(14) 



with A 



^) and B 



C J2 



where J„ = 



J„(Eo). Note that (fT^ correctly yields the equipartition 
for P oo : hHi = T J^. The expression (|12|) is more 
illuminating if we use a spectral representation of the 
Green function to transform the discrete sum over the 
Matsubara frequencies in an integral : 



hHi = [ —hujp{ijj)fB{uj) 

J-oo 



(15) 



p{uj) = -juj / lmGciq,ujn -ii^ + 0+) (16) 



tOJ 



(cg2--i^a;2 + Eo + /^(a;))2 + (/^'(u;))2 

where G';:^((j,w„) Z^b G'ab^(9, w„), p{uj) is the density 
of states and fsi^) the Bose factor. Eq. ((T5|l is simply 
the internal energy of free excitations, whose density of 
states is given by p(w), computed self-consistently within 
the variational method. As p{ijj) does not depend, at this 
order, on temperature, all the temperature dependence is 
contained in the Bose factor. We thus obtain the specific 
heat 



a = m? 



-t-oo 



du) p{lu)uj'^ 
47r sinh^ f3nuj/2 



+ 0{h) (17) 



Due to the Bose factor, the low T behavior is governed 
by the low lu behavior of the density of states. From 
(|14|l we see that p{uj) ^ ~ujIq{uj) ~ w^, which leads to 
Cv{T) ~ in all dimensions d > 2. Surprisingly, the 
linear term in the specific heat cancels at lowest order 
in h. This cancellation, and the resulting behavior 
occurs in fact to all orders in fi, as is shown below. This 
results in the following behavior at low temperature 



0{{T/hY) (18) 



Ci will be given in jSJj and we discuss here the dominant 
contribution (|17|l . Although the behavior is indepen- 
dent of the dimension 2 < d < 4, the coefhcient itself 
does depend on it. More precisely 



o{h,{i3hy 




in" 



Rr 



(19) 



where Kd — Sd/ (Stt)'', Sd being the volume of the sphere 
in dimension d and Fc^ (x) a scaling function. Choosing a 
spherical Brillouin zone, Fc^{x) — 2/7r(/2(a;)^//3(a;))^/^, 
with fnix) = J^^^ duu'^~^ I {v? -\- 1)". If we compare this 
result to the Debye law for pure phonons. 



aDcbyc ^Kd 



(20) 



we see that d = 3 appears as a particular dimension 
below which the specific heat is lowered by the disorder 
although it is enhanced above. In three dimensions it 
remains of the Debye form with a prefactor governed by 
the scaling function Fc^{x): 



Fc, (x) - 1 a; > 1 
FcA^) 5— a; x-^l 



(21) 
(22) 



Note that since one does recover exactly the Debye law at 
vanishing disorder, the corrections at weak disorder are 
rather small. At stronger disorder i?c ~ a the specific 
heat is significantly lowered by the disorder. 

In the case d < 2, it was noticed in TT) that to obtain 
the correct behavior of dynamical quantities (e.g. the 
conductivity cr(w) ~ w^) through linear response within 
the Matsubara equilibrium approach, one should enforce 
the MC {Tj), which arises naturally within the dynami- 
cal methods |Q. If we use this condition we obtain for 
Cy{T) the expression H17|l where in the numerator ■§iLlP' 
is replaced by — 2(i-d) ^0 ^^'^ ^^^^ leads to a neg- 

ative specific heat at low T. However, if we use the ther- 
modynamic condition, i.e. compute Uc by minimizing the 
variational free energy, the internal energy can again be 
written as (|12|l which guarantees the equipartition and a 
positive Cy{T). Then J(a;„) ~ w^, and Cv{T) vanishes 
exponentially at low T, a signature of a gapped excitation 
spectrum. This apparent discrepancy between thermo- 
dynamic and dynamic quantities within the GVM, which 
occurs only for d < 2, is a hint of possible complications 
due to slow equilibration in these systems, as observed 
e.g. in Coulomb glass systems f |3^l33j l. 

Extracting the T-dependence to all orders in fi is quite 
difficult, but we found a drastic simplification of the tem- 
perature dependence, at order T^, of t integrals such as 
in Q and (|5J). Schematically, for any polynomial form 
V{x) — xP, expressing B^ and B{tY as multiple Matsub- 
ara sums using (|10|l and performing the t integral, one 
explicitly checks that cancellations occur term by term 
in the multinomial expansion of B{tY ■ Details are given 
in '33|. This further allows to derive the cancellation of 
the contribution linear in temperature of to all orders. 



4 



arising naturally as a consequence of the marginality con- 
dition ((Tj), which, here for d> 2, holds automatically. We 
find that the self-energy takes the interesting form: 

S + /K) = S + A/(l - Sno) + /K) (23) 
E + M = + 0(T4) (24) 
= M + T^i^^^ i^n) + ■■■ (25) 

with S = + 0(r2) and /(w„) - |cj„| + 0(cj2). The 
absence of T^, dependences in Eq. (|24I25(I allows to 
show, using marginality, that no nor term survives 
in (O for d>2. For d < 2 the same conclusion as above 
persists (simply replacing Sq by E^"^). 

The drastic simplification in T dependence found here 
is also useful to compute analytically the putative linear 
term in C„ (T) in the quantum spin glasses. Both models 
in and 0| are solved via a one step RSB solution, 
hence analogous to our d <2 case, and thus are gapless 
provided the marginality condition is enforced (see dis- 
cussion above). Ref. studies the specific heat of the 
quantum p-spin model and obtain numerical evidence for 
a linear contribution. We find that this term vanishes. 
In Ref. the SU(N) fully connected Heisenberg spin 
glass (of spin S) was solved for infinite TV. There too 
it was found that the dominant term (analogous to Co 
here) in a \/S semiclassical expansion also behaves as 
~ . Whether it holds to higher order in 1/5 has been 
discussed in [s^ and [15| and is currently under investi- 
gation js^l- In all cases a structure similar to H23I24I25|) 
also holds. 

One can discuss the validity of the mean field method. 
The GVM can also be used to describe non periodic sys- 
tems (e.g . directed polymer) in terms of a power law 
V{B) II2.I3H and then becomes exact for m — > 00 with 
R{u) = -mV{u'^ /m). We thus expect that the be- 
havior is exact in this limit in the full RSB case. For the 
periodic problem it is always an approximation which 
describes the vibrational phonon-like excitations of the 
pinned system and may not treat accurately excitations 
such as solitons. Whether these solitons or other type of 
excitations such as dislocations for a lattice could reestab- 
lish a linear contribution to the specific heat from two 
level systems type behavior remains to be investigated. 

This weak T-dependence of the specific heat in the 
quantum regime is relevant to the so-called vortons 
modes of the vortex lattice of superconductors. Other 
contributions such as core electrons j^, or their inter- 
action with vortons (responsible for the dissipation in 
the quantum dynamics of the vortons) yield linear spe- 
cific heat |26| which may dominate over the present one 
in the quantum regime. Whether this also holds in the 
classical regime for the vortons is examined in ^] . 

In conclusion we have obtained from first principles the 
specific heat of an elastic quantum pinned system, such 
as the Bragg glass phase of a lattice in presence of impu- 
rities. We found that ^ in c? = 2, 3 due to a non 



trivial mechanism of cancellation of the linear term. This 
simplification which occurs from the marginality condi- 
tion will be useful for a general understanding of the spe- 
cific heat behavior within the mean field theory. The 
question remains concerning the contribution of other 
types of excitations such as plastic deformations. Fur- 
ther analytical and numerical investigations would be de- 
sirable to a better understanding of the behavior of the 
specific heat in glasses. 
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